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Sea  Surface  Multipath  Effects  on  Ship  Radar 
Radiated  Power  Determination 


L  INTRODUCTION 

One  of  the  factors  in  defining  a  ship's  radar  is  effective  radiated  power 
(ERP).  Knowing  the  ERP  of  a  radar  along  with  scan  rate,  frequency,  and  beam 
width  characterizes  the  radar  in  terms  of  type,  size,  and  transmitter  power. 

This  information  can  be  useful  in  quantifying  ship  traffic  in  port  areas. 

To  obtain  accurate  ERP  data  from  ship  radars  over  a  short  5-20  nautical 
mile  link,  multipath  effects  must  be  accounted  for.  Multipath  effects  have  been 
studied  extensively  [1-12].  Sea  surface  multipath  results  from  two  different 
phenomena  -  optical  interference  and  ducting.  Optical  interference  occurs  when 
the  radar  signal  from  the  ship  to  the  station  travels  in  both  a  direct  path  and  a 
path  reflecting  off  the  water  (Fig.  1).  The  signals  from  the  two  paths  combine  in 
and  out  of  phase  depending  on  the  distances  and  radar  frequency  and  cause  the 
signal  strength  to  vary.  The  reflected  signal  can  also  be  affected  by  sea  state 
and  the  presence  of  ship  wakes  and  other  sea  obstructions  (buoys  and  other 
ships). 


Fig.  1  Optical  Interference 


Multipath  from  optical  interference  can  be  specular  or  diffuse.  The 
specular  component  is  coherent  with  respect  to  the  direct  signal,  and  is  well- 
defined  in  terms  of  amplitude,  phase,  and  incident  direction.  The  diffuse 
component  has  a  random  nature,  and  arises  from  scattering  sources  from  many 
directions  [13].  Generally,  specular  multipath  is  modeled  with  the  diffuse 
component  added  in  through  a  reflection  coefficient.  The  diffuse  component 
becomes  stronger  at  low  grazing  angles  and  for  high  sea  states.  With  diffuse 
scattering,  the  deep  nulls  associated  with  optical  interference  effects  are 
smoothed  out. 

The  second  multipath  phenomenon  is  ducting.  A  duct  is  a  layer  of  air  with 
different  properties  than  its  surroundings  which  acts  as  a  waveguide  to  trap 
electromagnetic  energy.  These  refractive  gradients  in  the  atmosphere  can  cause 
over-the-horizon  fields  to  be  tens  of  decibels  higher  than  expected.  Refractivity 
profiles  can  be  measured  by  radiosondes  to  indicate  the  presence  of  ducts.  This 
technique  is  not  suitable  in  many  situations  since  the  radiosonde  emits  high 
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power  and  must  be  flown  at  a  range  of  heights.  Also,  the  radiosonde  only 
measures  the  refractivity  profile  directly  above  the  receiving  station,  whereas, 
depending  on  the  transmitter  location,  it  may  be  the  profile  several  kilometers 
away  which  is  affecting  the  signal.  Refractivity  profiles  change  on  a  hourly 
basis,  so  obtaining  adequate  refractivity  data  to  predict  ducting  can  be 
unfeasible  for  other  than  short  term  experiments. 

Ducting  can  occur  with  various  meteorological  conditions.  Evaporation 
ducts,  caused  by  the  rapid  decrease  in  humidity  with  increasing  altitude  just 
above  the  ocean's  surface,  occur  with  the  base  of  the  duct  at  the  sea  surface. 
Surface-based  ducts  generally  occur  at  higher  heights  than  evaporation  ducts. 
They  are  caused  by  temperature  and  humidity  inversions  when  a  warm  dry  air 
mass  moves  over  a  cool,  moist  air  mass.  [15]  Ducting  effects  are  more 
prominent  at  longer  ranges  and  allow  reception  of  RF  signals  at  much  longer 
ranges  than  would  normally  occur.  Fig.  2  shows  a  schematic  of  ducting. 
Appendix  A  gives  more  detailed  information  about  the  ducting  phenomena. 
Appendix  B  summarizes  common  characteristics  of  navigation  radars. 


Fig.  2  Ducting 

Various  methods  can  be  used  to  counteract  the  effects  of  multipath.  For 
communication  links,  the  use  of  spatial  and  frequency  diversity  has  been  used. 
Kuhnert  and  Gelerman  [14]  describe  the  use  of  two  frequencies  and  two  heights 
on  a  receiving  tower  to  counteract  ducting.  An  adaptive  equalizer  was  used. 

This  study  is  broken  into  six  sections.  The  first  describes  various 
computer  codes  used  for  analyzing  multipath.  The  second  describes  details  of 
the  ERP  determination  process.  In  the  third  section,  a  fixed  sea-surface  link  is 
used  to  analyze  sea-surface  multipath,  including  tidal  effects.  The  fourth 
section  describes  upper  air  soundings  performed  over  the  Chesapeake  Bay.  The 
sounding  data  provides  refractivity  profiles  that  are  used  as  input  to  a 
computer  code  to  further  analyze  multipath  and  particularly  ducting  effects  on 
ERP  data  collection.  The  fifth  section  describes  models  to  predict  radio 
frequency  (rf)  attenuation  due  to  rain.  Along  with  multipath  effects,  rain 
attenuation  can  also  affect  accurate  ERP  determination.  The  last  section  gives 
results  of  the  ERP  data  collection  experiment  which  involved  taking  data  on 
several  hundred  ships.  Several  techniques  are  used  to  estimate  the  accuracy  of 
the  data. 
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II.  COMPUTER  CODES  FOR  MULTIPATH  MODELING 

Several  computer  codes  are  available  to  model  sea  surface  multipath 
effects.  The  EREPS  (Engineer's  Refractive  Effects  Prediction  System)  code  [15,16] 
models  optical  interference,  diffraction,  tropospheric  scatter,  refraction, 
evaporation  and  surface-based  ducting,  and  water  vapor  absorption  under 
horizontally  homogeneous  atmospheric  conditions.  The  refractivity  profile  of 
the  atmosphere  is  approximated  by  the  use  of  an  effective  earth  radius 
(Appendix  C).  Wave  height  is  modeled  as  a  function  of  wind  speed  as  described 
in  Appendix  D.  Figure  3  shows  a  sample  output  of  the  EREPS  code,  showing  a 
comparison  of  propagation  loss  versus  range  with  and  without  a  13  m  high 
evaporation  duct.  In  this  case,  the  presence  of  a  duct  greatly  decreases  the 
propagation  loss  for  distances  beyond  20  nm. 


Rang©,  nautical  miles 

Fig.  3  Propagation  loss  versus  range  for  radar  at  9400  MHz  with  and  without 
evaporation  duct  (EREPS  code,  Wind  speed  =  10  knots,  Transmitter  height  =  60 
ft,  Receiver  height  =180  ft,  no  ducting,  omni  transmitter) 

Several  EREPS  runs  comparing  results  with  varying  wind  speed,  ducting 
conditions,  receiver/transmitter  heights,  and  frequencies  are  presented  in 
Appendix  E.  Runs  were  also  done  using  spatial  diversity  (where  signals  from 
two  receiver  heights  are  combined).  The  multipath  lobing  dramatically 
decreases  for  this  case.  It  was  found  that  the  EREPS  code  is  only  good  for  coarse 
predictions  because  the  refractivity  profile  is  assumed  to  be  linear.  The  duct 
features  are  approximate  as  well. 
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The  Northam  Multipath  code  [17,18]  models  multipath  effects  for 
microwave  radar  signals  for  low  grazing  angle,  forward  scatter,  over-water 
multipath.  Northam  characterizes  over-water  multipath  as  having  both 
specular  and  diffuse  components.  These  components  are  influenced  by  sea 
state,  reflectivity  of  the  water,  radar  frequency,  and  polarization.  In  Northam's 
model,  the  sea-state  is  represented  by  the  root  mean  square  (RMS)  surface 
displacement  above  the  mean.  The  specular  component  is  treated  as  a 
deterministic  component,  with  parameters  unknown  to  the  tracking  system. 

The  diffuse  component  is  treated  as  stochastic.  This  code  was  specifically 
designed  to  model  multipath  effects  on  pulse  structure.  It  is  a  Fortran  code 
which  currently  runs  on  a  Macintosh  computer. 

The  Radio  Physical  Optics  (RPO)  [19]  code  is  similar  to  EREPS  but  includes 
the  use  of  refractivity  profiles.  It  can  calculate  range-dependent 
electromagnetic  (EM)  system  propagation  loss  for  a  heterogeneous  atmosphere. 
The  radio-frequency  index-of-refraction  can  be  varied  both  vertically  and 
horizontally.  This  code  was  developed  by  the  Propagation  Division  at  NOSC  in 
San  Diego,  California,  which  has  developed  several  other  propagation  codes  and 
releases  them  over  the  Internet.  The  RPO  code  was  found  to  be  useful  for  much 
of  the  modeling  required  in  this  report.  Its  main  limitation  was  in  needing  large 
amounts  of  refractivity  data,  which  change  on  an  hourly  basis  and  are  difficult 
to  obtain.  In  our  case,  using  one  refractivity  profile  for  the  entire  9  nm  over¬ 
water  link  did  not  provide  exact  predictions  for  link  performance,  although 
relative  performance  estimates  could  be  obtained  and  are  valuable  for 
analyzing  the  link  data.  Walters,  et.  al.,  [20]  has  used  the  RPO  code  for 
validating  radar  performance  using  refractivity  data  from  a  tethered  radiosonde 
with  some  success. 

The  Advanced  Propagation  Model  (APM  1.0)  from  NOSC  is  a  hybrid  model 
that  consists  of  four  sub-models:  flat  earth,  ray  optics,  extended  optics, 
and  split-step  parabolic  equation  (PE).  APM  effectively  merges  both  the 
Radio  Physical  Optics  (RPO)  model  and  the  Terrain  Parabolic  Equation 
Model  (TPEM).  The  result  is  a  new  and  improved  EM  propagation  model  that 
can  model  range-dependent  refractivity  environments,  variable  terrain,  range- 
varying  dielectric  ground  constants  for  finite  conductivity  and  vertical 
polarization  calculations,  troposcatter,  and  gaseous  absorption.  This  code  has 
only  recently  become  available  for  public  distribution  and  was  not  used  in  this 
study. 


A  listing  of  other  computer  codes  for  multipath  modeling  is  given  in 
Appendix  F.  Results  of  computer  modeling  will  be  compared  to  experimental 
data  in  the  following  sections. 
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III.  ERP  Calculations 


The  ERP  data  collected  in  this  report  was  taken  at  the  Chesapeake  Bay 
Detachment  (CBD)  of  NRL  at  Chesapeake  Beach,  Maryland  [21].  Collection 
antennas  are  set  118  ft.  above  the  water  level  on  a  cliff  overlooking  the  bay. 

The  x-band  (8-12  GHz)  and  s-band  (2-4  GHz)  antennas  are  omnidirectional 
bicones  that  are  shielded  in  the  rear  by  an  enclosure,  giving  them  a  pattern  that 
is  relatively  flat  to  +/-  60°.  This  type  of  pattern  enables  the  signal  amplitudes 
to  be  calculated  with  just  one  value  for  the  antenna  gain,  while  rejecting  signals 
in  other  directions.  A  receiver  system  attached  to  the  antennas  allowed  range, 
frequency,  scan  rate,  and  illumination  time  (pulse  width)  data  to  be  taken  along 
with  signal  amplitude. 

The  space  path  loss,  L,  of  power  of  a  radio  wave  in  space  is  [22]: 


f4jtR) 

2 

facfR'l 

v  x  ) 

V  C  J 

where 


R  =  range 
X=  wavelength. 

For  radar  signals  from  ships,  optical  interference  from  the  sea  surface 
creates  a  multipath  effect,  where  the  signal  strength  varies  as  the  range  is 
varied.  Figure  1  shows  the  case  where  a  signal  from  a  ship  travels  directly  to 
the  receiver  and  also  takes  a  path  which  bounces  off  the  water.  The  signals 
from  the  two  paths  combine  in  and  out  of  phase  depending  on  the  distances  and 
frequency,  so  the  signal  strength  will  vary  as  the  ship  moves,  as  the  tide  varies, 
and  as  the  sea  state  changes.  This  phenomenon  is  also  affected  by  the  presence 
of  ship  wakes  and  other  sea  obstructions  (buoys  and  other  ships).  Figure  4 
shows  an  example  of  interference  effects  for  an  over-water  link  and  two 
receiver  heights  using  the  EREPS  code.  As  the  receiver  (or  transmitter)  height 
changes,  the  interference  pattern  shifts  so  that  the  peaks  and  nulls  occur  at 
shifted  locations. 
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Range,  nautical  miles 

Fig.  4  Propagation  Loss  versus  range  for  over-water  link  at  9400  MHz  with 
varying  receiver  heights  (EREPS  code,  Wind  speed  =  10  knots,  Transmitter 
height  =  60  ft,  no  ducting,  omni  transmitter) 


Figure  5  shows  an  example  of  wind  speed  effects  for  an  over-water  link 
using  the  EREPS  code.  Wave  height  is  modeled  as  a  function  of  wind  speed  in 
the  EREPS  code  (Appendix  D).  High  winds  correlate  to  large  wave  heights,  which 
mute  out  the  interference  pattern  and  reduce  the  null  depth.  The  peaks  of  the 
signal  are  also  lowered.  At  5  nm,  the  peak  of  the  signal  is  lowered  by  1  dB  with 
a  15  knot  wind  (2.5  foot  wave  height). 
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Range,  nautical  miles 

Fig.  5  Propagation  Loss  versus  range  for  Choke  Point  Monitoring  System  at 
9400  MHz  at  two  wind  speeds  (EREPS  code,  Transmitter  height  =  60  ft,  Receiver 
height  =180  ft,  no  ducting,  omni  transmitter) 

Ship  blockage  and  multipath  off  ship  structures  can  have  a  considerable 
effect  on  signal  reception.  This  will  be  direction  dependent.  Rocking  of  a  ship 
due  to  high  waves  and  swells  causes  radar  signals,  with  their  narrow  beams,  to 
fade  in  and  out.  In  this  case,  the  direction  of  the  waves  (usually  correlated  with 
wind  direction)  with  respect  to  the  ship  and  the  receiver  is  a  big  factor. 

For  x-  and  s-band,  optical  interference  effects  dominate  out  to  the  radio 
horizon  for  ground-to-ground  transmission.  For  an  x-band  radar  with  a  100 
foot  transmitter  and  receiver,  this  is  ~  18  nm.  This  area  is  called  the  optical 
region.  After  that,  ducting  and  refractive  index  phenomena  will  have  the 
dominant  effect.  For  ground- to-space  transmission,  interference  effects  will  be 

noticeable  at  elevation  angles  higher  than  ~  1°. 

For  the  ERP  determination,  the  radar  signal  amplitude  was  plotted  against 
range  along  with  a  free  space  loss  curve.  The  peaks  of  the  signal  should  be  6  dB 
above  the  free  space  curve  for  the  case  where  no  waves  are  present.  This 
decreases  to  ~5  dB  for  a  more  realistic  case  where  wave  heights  are  between 
1.5-3  feet.  For  the  calculations  in  this  report,  5  dB  was  used.  The  offset  needed 
to  place  the  free  space  curve  5  dB  below  the  peaks  of  the  data  was  the  ERP  of 
the  signal.  Figures  6  and  7  show  example  plots  for  s-band  and  x-band  radars. 
Other  data  plots  are  presented  in  Appendix  G. 
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Using  the  computer  models,  it  was  found  that  for  a  receiver  height  of 
116  feet  refractivity  deviations  from  a  standard  atmosphere  affect  ERP 
predictions  for  data  beyond  approximately  8  nm.  This  is  because  at  longer 
ranges,  the  angle  the  signal  travels  through  the  atmosphere  is  smaller  and  so 
atmospheric  bending  is  more  pronounced.  Strong  duct  conditions  can  affect  ERP 
predictions  to  within  4-5  nm.  Our  ERP  predictions  are  normally  taken  from  data 
4-6  nm  in  range,  so  should  only  be  affected  by  strong  ducting  conditions  that 
occur  mainly  in  the  afternoons  of  summer  months  with  still  air. 

The  results  of  the  EREPs  model  can  be  used  to  determine  the  transmitter 
height.  The  frequency,  receiver  height  and  wave  heights  are  inputted  to  the 
code  and  the  transmitter  height  varied  until  the  multipath  lobes  align. 


Fig.  6  S-band,  3053  MHz,  Radid  328179,  in-bound,  11/24/97,  6:20  am, 
Transmitter  height  115  feet,  receiver  height,  116  ft.  Free  space  loss  is 
approximately  5  dB  below  the  peak  signal  (varies  slightly  with  sea  state).  The 
offset  for  the  free  space  loss  corresponds  to  the  ERP  of  the  generating  signal,  in 
this  case  104  dBm. 
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Fig.  7  X-band,  9384  MHz,  Radid  142763,  out-bound,  5/22/97,  1200, 

Transmitter  height  108  feet,  receiver  height,  116  ft.  ERP  =  105  dBm. 

Once  the  ERP  is  known,  transmitter  power  can  be  calculated  by  knowing 
the  vertical  and  horizontal  beamwidths.  The  vertical  beamwidth  (vbw)  is 
estimated  to  be  ~  25°.  The  horizontal  beamwidth  (hbw)  is  determined  by 
dividing  the  illumination  time  by  the  scan  time  and  multiplying  by  360°.  The 
antenna  gain,  G,  in  dB  (including  efficiency,  r|)  is  then 

G=10LOG( — — )  for  hbw  and  vbw  in  radians. 

hbw  *  vbw 

27000 

G  ~10  LOG( - )  for  hbw  and  vbw  in  degrees. 

hbw  *  vbw 

10* 

The  transmitter  power,  P,  in  kW  is 

where  x  =  ((ERP-G+L)/10)  and  L  is  the  line  loss  from  the  transmitter  to  the 
antenna,  approximately  0.8  dB. 
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IV.  STATIC  MULTIPATH  EFFECTS 


A  stationary  radar  situated  on  a  100  foot  tower  on  Tilghman  Island  was 
used  as  a  benchmark  for  analyzing  multipath  effects  across  the  Chesapeake  Bay 
from  CBD.  The  distance  across  the  Bay  at  this  site  was  ~  10  nm.  Data  taken  on 
this  radar  could  be  correlated  with  tidal  data  and  refractivity  data  from  upper 
air  soundings.  Much  of  the  signal  variations  could  be  attributed  to  changes  in 
tide  height.  However,  particularly  during  the  summer  months,  there  were  large 
signal  variations  that  were  primarily  due  to  ducting.  With  the  tidal  variations, 
the  signal  level  followed  approximately  an  11  hour  period,  whereas  with  the 
ducting  variations,  the  signal  fluctuations  occurred  with  periods  of  a  few 
minutes.  Figs.  8  &  9  shows  the  signal  following  tidal  variations  over  a  period  of 
two  days  and  ten  days.  Tidal  data  was  obtained  from  the  Next  Generation 
Water  Level  Measurement  System  (NGWLMS),  produced  by  the  National  Ocean 
Service  (NOS)  (NOAA).  The  data  (height  above  Mean  Lower  Low  Water)  was 
taken  at  Naval  Academy  in  Annapolis,  20  nm  north  of  CBD,  but  adjusted  ahead 
in  time  by  1  hour  14  minutes  to  match  the  tide  levels  at  CBD  (based  on  tide 
charts).  Mean  Lower  Low  Water  is  a  tide  reference  point  relative  to  survey 
markers  at  the  Anapolis  site.  Fig.  10  shows  tidal  variations  over  a  period  of  30 
days. 
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Fig ♦  8  Tilghman  Island  radar  amplitude  versus  tide  variations  in  the 

Chesapeake  Bay 
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Time  (days) 

Fig.  10  Tidal  variations  beginning  8/23/95  at  Naval  Academy  in  Annapolis 
(from  the  Next  Generation  Water  Level  Measurement  System  (NGWLMS),  from 
the  National  Ocean  Service  (NOS)  (NOAA))  MLLW  (Mean  Lower  Low  Water) 


Some  fluctuations  could  be  attributed  to  ship  wakes.  The  signal  level 
variations  due  to  tidal  changes  (optical  interference  effect)  could  be  modeled 
with  the  EREPS  code,  however,  it  was  found  that  the  refractive  index  could 
change  significantly  from  day  to  day  and  no  one  single  value  could  be  used  to  fit 
all  the  data.  Fig.  1 1  shows  the  predictions  of  propagation  loss  versus  tidal 
height  using  the  EREPS  code  using  a  standard  refractive  index.  Fig.  12  is  similar 
except  the  calculation  uses  higher  wind  speed  (which  corresponds  to  higher 
wave  height).  As  discussed  earlier,  increased  wave  height  decreases  the 
severity  of  the  multipath  variations. 
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Propagation  Loss  (dB) 


Fig.  11  Propagation  Loss  versus  relative  tidal  height  for  Choke  Point 
Monitoring  System  at  9404.5  MHz,  EREPS  prediction  for  Tilghman  Island  radar 
at  various  heights  (Wind  speed  =  10  knots,  no  ducting,  omni  transmitter) 
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Propagation  Loss  (dB) 


Fig.  12  Propagation  Loss  versus  relative  tidal  height  for  Choke  Point 
Monitoring  System  at  9404.5  MHz,  EREPS  prediction  for  Tilghman  Island  radar 
at  various  heights  (Wind  speed  =  20  knots,  no  ducting,  omni  transmitter) 
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V.  UPPER  AIR  SOUNDINGS  AND  CORRELATION  TO  MODELING 


Upper  air  soundings  on  the  Chesapeake  Bay  between  Tilghman  Island  and 
CBD  were  performed  on  three  days  using  a  tethered  aerostat  and  a  remote 
temperature/pressure/humidity  module  (Fig.  13).  Two  to  three  refractivity 
profiles  were  obtained  per  day.  The  data  was  used  with  the  RPO  and  EREPS 
codes  to  correlate  with  radar  data.  It  was  also  used  with  the  RPO  code  to  help 
validate  the  ERP  data. 


Fig.  13  Tethered  aerostat  (a)  and  remote  sensing  module  (b)  used  for  upper- 
air  soundings 

For  the  experiment,  the  sensing  module  (ENV-50-HUM  from 
Sensormetrics)  was  connected  below  the  aerostat  in  a  well-ventilated  shielded 
container.  It  was  necessary  to  keep  direct  sunlight  off  the  temperature  sensor. 
The  aerostat  was  reeled  out  on  a  marked  rope  in  increments  of  10  feet  every 
minute.  The  time  delay  was  needed  because  the  sensing  module  was  designed 
to  take  one  measurement  per  minute,  in  part  due  to  the  settling  time  of  the 
sensors.  The  height  increments  were  increased  to  20  feet  after  ~  150  feet,  to 
50  feet  after  500  feet,  and  to  100  feet  after  1500  feet.  The  aerostat  lift  only 
allowed  the  sensor  to  reach  2000  feet.  One  set  of  measurements  took  about  one 
hour  to  perform.  The  aerostat  was  then  reeled  in  and  the  temperature, 
pressure,  and  humidity  data  extracted  from  the  memory  of  the  sensing  module 
with  a  PC  interface. 
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Equations  for  converting  temperature,  pressure,  and  humidity  data  into 
refractivity  and  modified  refractivity  are  given  in  Appendix  C.  Barometric 
pressure  was  used  to  determine  height  (see  Appendix  C),  however,  rope 
measurements  were  also  used  at  the  lower  heights  to  get  more  accurate 
readings. 

Fig.  14  shows  modified  refractivity  profiles  taken  in  the  middle  of  the 
Chesapeake  Bay  during  a  warm  day  with  light  winds  (June  17,  1996).  Modified 
refractivity  is  plotted  because  it  shows  the  effects  of  ducts  more  clearly  than 
when  using  refractivity.  No  strong  surface  ducts  were  present,  though  there 
was  an  evaporation  duct  present  at  the  1:30  pm  sounding.  Fig.  15  shows  the 
RPO  modeling  done  using  a  standard  atmosphere  refractivity  profile.  Fig.  16 
shows  the  RPO  modeling  done  using  the  10:30  am  modified  refractivity  profile. 
As  can  be  seen,  the  interference  pattern  structure  (used  to  derive  ERP 
measurements)  starts  breaking  up  after  ~  10  km  for  transmitter  and  receiver 
heights  of  30  m. 


Modified  Refractivity,  M-units 

Fig.  14  Chesapeake  Bay,  June  17,  1996.  Warm  day  with  light  winds,  no 

strong  ducts  present. 
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Fig.  15  RPO  modeling  using  refractivity  data  from  stdatm,  0  350,  1000.  468 


110  115  120  125  130  135  140  145  150 
Propagation  Loss,  dB 

Fig.  16  RPO  modeling  using  refractivity  data  from  June  17,  1996  Set  #1, 
10:30  am,  with  a  tide  height  of  1  ft.  and  transmitter  height  of  30  m. 
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Fig.  17  shows  modified  refractivity  profiles  taken  during  a  hot,  calm  day 
on  the  Chesapeake  Bay.  For  this  day,  strong  evaporation  and  surface-based 
ducts  were  present.  Fig.  18  shows  the  RPO  modeling  results  using  the  12:20  pm 
modified  refractivity  profile.  In  this  case,  the  interference  pattern  structure 
starts  breaking  up  after  ~  8  km  for  transmitter  and  receiver  heights  of  30  m. 


Modified  Refractivity,  M-unite 

Fig.  17  Chesapeake  Bay,  July  16,  1996.  Middle  of  hot,  calm  day  shows  strong 

elevation  and  surface-based  ducts. 

Radiosondes  were  also  considered  for  performing  these  measurements. 

The  cost  for  procuring  a  radiosonde  system  turned  out  to  be  prohibitive.  Also, 
the  radiosonde  would  still  need  to  be  tethered  to  obtain  accurate  data  readings 
at  low  altitudes  to  measure  evaporation  and  surface-based  ducts  because  the 
sensor  reading  needs  time  to  stabilize.  Normally  a  radiosonde  is  released  with  a 
balloon  and  measurements  taken  every  100-200  feet.  The  radiosonde  can 
measure  atmospheric  conditions  at  heights  up  to  80,000  feet.  Radiosonde 
measurements  are  taken  on  a  regular  basis  around  the  country,  primarily  near 
airports,  as  an  aid  to  navigation.  Much  of  the  radiosonde  data  is  available  over 
the  internet. 
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Fig.  18  RPO  modeling  using  refractivity  data  from  July  16,  1996  Set  #2, 

12:20  pm 


Fig.  19  shows  modified  refractivity  profiles  taken  during  a  cool  windy  day 
on  the  Chesapeake  Bay  (Sept.  19,  1996).  Fig.  20  shows  the  RPO  modeling  done 
using  this  modified  refractivity  profile.  The  temperatures  were  mild,  near  70°F, 
with  12-15  knot  winds,  1.5-2  foot  waves,  and  clear  skies.  For  this  day,  no 
strong  ducts  were  present.  This  is  to  be  expected  because  of  the  weather 
conditions.  The  cooler  air  will  not  promote  as  much  evaporation  of  water  (to 
create  evaporation  ducts),  and  the  high  winds  mix  up  the  air,  preventing  any 
strong  surface-based  ducts  from  forming.  The  refractivity  patterns  remained 
similar  for  the  three  data  sets.  The  three  curves  shift  over  time.  This  is 
because  of  the  rise  in  temperature  in  the  afternoon,  which  corresponds  to  a  fall 
in  relative  humidity.  The  curves  still  maintained  the  same  slope,  though,  which 
indicates  the  effective  earth  radius  value  remained  constant  (Appendix  C),  and 
the  propagation  through  the  atmosphere  should  not  change  significantly. 

Temperature  and  relative  humidity  graphs  of  the  soundings  along  with 
other  data  are  presented  in  Appendix  H. 
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Modified  Refractivity,  M 


Fig.  19  Chesapeake  Bay,  Sept.  19,  1996.  Sunny,  mild,  windy  day 
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110  115  120  125  130  135  140  145  150 
Propagation  Loss,  dB 

Fig.  20  RPO  modeling  using  refractivity  data  from  Sept.  19,  1996  Set  #2, 

12:00  am 

Figure  21  shows  the  signal  amplitudes  from  the  Tilghman  Island  radar 
during  the  same  period  as  when  the  July  16,  1996  soundings  data  was  obtained. 
As  can  be  seen,  signal  levels  fluctuated  by  +/-  5  dB  during  this  period, 
uncorrelated  to  tide  height.  Attempts  to  make  a  direct  prediction  of  signal 
levels  to  the  RPO  code  failed.  As  shown  in  Fig.  17,  the  refractivity  profile  can 
change  dramatically  in  1.5  hours.  Along  with  observed  changes  in  the 
refractivity  profile  on  the  order  of  minutes,  the  assumption  that  the  profile  will 
be  constant  across  the  entire  rf  path  (across  the  Bay)  is  also  incorrect,  and  this 
contributes  to  the  inability  to  use  the  refractivity  data  and  the  RPO  code  to 
make  accurate  predictions.  However,  the  RPO  code  can  give  insight  into  why 
such  dramatic  signal  level  variations  are  seen  in  the  Tilghman  radar  data. 

When  small  changes  are  input  to  the  refractivity  data  in  the  RPO  code,  the  nulls 
and  peaks  of  the  signal  level  pattern  move.  This  movement  is  more  pronounced 
at  longer  distances.  Using  the  technique  of  varying  refractivity  data  in  known 
increments,  the  RPO  code  can  be  used  to  predict  the  magnitude  and  to  a  smaller 
extent  the  time  scale  of  the  signal  variation.  It  also  can  predict  the  distance  and 
transmitter  height  levels  where  the  standard  interference  pattern  of  the  signal 
levels  breaks  down  due  to  various  refractivity  conditions. 
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The  use  of  GPS  or  remote  sensing  satellites  to  predict  refractivity  profiles 
over  large  areas  may  provide  a  feasible  method  of  obtaining  refractivity 
profiles  over  links  in  the  future,  so  that  prediction  methods  such  as  the  RPO 
code  could  be  used  more  fruitfully.  The  GPS  technique,  called  the  direct 
inference  technique  (DIT),  compares  the  observed  interference  pattern  as  the 
satellite  moves  through  low  elevation  angles  with  patterns  predicted  with 
known  refractivity  profiles.  The  profile  with  the  best  fit  is  then  assumed  to  be 
correct.  [23-25] 
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Fig.  21  Average  Signal  Level  from  Tilghman  Island  radar  compared  to  tide 
height  for  July  16,  1996.  (note:  Time  in  days  in  zulu  time.  Time  for  sets  is  EST 
(Eastern  Standard  Time). 


Figure  22  shows  the  signal  amplitudes  from  the  Tilghman  Island  radar 
during  the  same  period  as  when  the  Sept.  19,  1996  soundings  data  was 
obtained.  As  can  be  seen,  signal  levels  fluctuated  by  only  +/-  2  dB  during  this 
period,  and  levels  can  be  partially  correlated  to  tide  height.  The  signal  levels 
are  much  more  stable  than  for  the  July  16,  1996  case,  which  is  expected,  since 
the  refractivity  profiles  for  the  September  case  do  not  show  significant  ducting 
or  variability  over  time. 


22 


H 

a. 

CD 

I 

CD 

c5‘ 


p 


Fig.  22  Average  Signal  Level  from  Tilghman  Island  radar  compared  to  tide 
height  for  Sept.  19,  1996.  (note:  Time  in  days  in  zulu  time.  Time  for  sets  is  EST 
(Eastern  Standard  Time). 


A  modeling  study  was  also  conducted  to  determine  the  height  where  the 
atmosphere  still  had  an  effect  on  a  20  nm  (or  less)  link  (with  transmitter  height 
of  30  m,  and  receiver  height  up  to  100  m).  The  slope  of  the  refractivity  profile 
was  drastically  altered  above  various  heights  and  any  changes  in  the 
propagation  loss  diagram  noted.  It  was  determined  from  this  study  that 
atmospheric  conditions  above  ~  500  m  had  almost  no  effect  on  the  20  nm  link. 
Atmospheric  conditions  from  50  to  500  m  had  a  minor  effect.  The  major  effect 
occurred  with  a  change  to  atmospheric  conditions  below  50  m.  From  the 
soundings  data  it  was  observed  that  evaporation  ducts  (which  occur  below  100 
m)  on  the  Bay  were  prevalent  during  the  two  summer  days  when  data  was 
taken.  The  size  and  shape  of  the  ducts  varied  greatly  in  the  space  of  a  few 
hours. 
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VI.  RAIN  ATTENUATION  MODELING 


Rain  attenuation  is  the  lessening  of  intensity  of  a  signal1  due  to  absorption 
and  scattering  by  raindrops  in  the  path  of  the  radar.  The  raindrops  make  it 
difficult  to  accurately  calculate  attenuation,  because  in  storms  the  size  and 
shape  of  raindrops  is  not  uniform  [26,27].  In  general,  attenuation  is  calculated 
based  on  average  rainfalls,  not  on  a  drop  to  drop  basis.  Rain  attenuation  is 
much  more  severe  for  x-band  than  s-band  radars,  in  part  because  the 
wavelength  of  x-band  (1"-1.5")  is  closer  in  size  to  the  diameter  of  a  raindrop. 

Crane’s  model  [28-31]  provides  a  simple  technique  for  estimating  rain 
attenuation  based  on  frequency,  rain  rate,  and  ground  distance.  Figs.  23  and  24 
graph  attenuation  versus  one-way  distance  for  x-band  and  s-band  radars. 

Other  rain  predictions  models  are  also  available  [32-34].  Appendix  I 
compiles  a  list  of  weather  data  sources  from  which  archived  and  real  time 
weather  data  can  be  obtained. 


Fig.  23  Crane  rain  model  prediction  of  x-band  signal  attenuation  versus  one¬ 
way  distance  for  various  rainfall  conditions 
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Fig.  24  Crane  rain  model  prediction  of  s-band  signal  attenuation  versus  one¬ 
way  distance  for  various  rainfall  conditions 

Attenuation  is  not  only  caused  by  rainfall.  Attenuation  occurs  when  any 
medium  interferes  with  the  radar’s  path  of  propagation.  This  includes  such 
things  as  atmospheric  gases,  clouds,  hail,  snow,  water  vapor,  fog  and  ice. 
However,  all  these  factors  have  less  attenuation  than  rain.  At  wavelengths 
greater  than  a  few  centimeters,  absorption  by  atmospheric  gases  is  generally 
thought  to  be  negligibly  small  except  where  very  long  distances  are  concerned 
[26].  Also  the  liquid  form  of  water  (rain),  causes  more  attenuation  than  its  gas 
(vapor)  or  solid  (snow,  ice,  hail)  form.  The  attenuation  caused  by  hail  is  one- 
hundredth  that  caused  by  rain,  ice  crystal  clouds  cause  no  sensible  attenuation, 
and  snow  produces  very  small  attenuation  even  at  the  excessive  rate  of  fall  of  5 
inches  an  hour  [27].  Clouds,  including  fog,  have  the  second  greatest  effect  on 
attenuation.  Fog  is  just  a  cloud  at  a  low  lying  level  in  the  atmosphere. 
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VII.  ERP  DATA  COLLECTION  AND  VALIDATION 


ERP  data  on  approximately  600  ships  traveling  on  the  Chesapeake  Bay 
was  obtained  over  a  period  from  June  1996  through  December  1997.  For  each 
individual  ship  radar,  power  level  data  and  range  was  obtained  from  the  CBD 
collection  system.  This  data  was  then  fed  into  a  computer  program  which  would 
plot  out  the  results  and  give  an  ERP  reading  based  on  the  closest  peak  on  the 
signal.  An  operator  would  then  examine  the  data  and  make  a  determination  as 
to  whether  an  accurate  ERP  could  be  obtained.  Occasionally  the  data  would 
contain  spurious  points,  very  few  points,  or  an  irregular  shape.  This  data  would 
be  discarded. 

Figures  25-28  show  histograms  of  ERP  level  and  frequency  versus  number 
of  radars  for  all  the  data  collected.  For  these  four  charts,  ships  with  multiple 
visits  were  counted  more  than  once.  ERP  levels  below  90  dB  represented  very 
marginal  data. 
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Frequency,  GHz 

Fig.  26  X-band  Radar  Freq.  Variation  from  CBD  data  (6/96-12/97) 
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Fig.  27  S-band  Radar  ERP  Variation  from  CBD  data  (June  1996-Dec.  1997) 
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Fig.  28  S-band  Radar  Frequency  Variation  from  CBD  data 
(June  1996-Dec.  1997) 

Two  methods  were  used  for  ERP  data  validation.  The  first  method 
involves  taking  data  from  a  known  radar  and  comparing  the  computed  ERP  with 
the  known  ERP.  So  far,  we  have  correlated  data  from  three  known  radars  and 
obtained  good  agreement  on  signal  levels.  Table  I  shows  the  results 
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Table  I.  Validation  Points  for  ERP  Level  Data 


Radid 

Radar  Type 

ERP  Level,  dB 
CBD/Listing 

Power,  kW 
CBD/Listing 

185489 

Furuno  8505DA 

97. /98. 

4. 4/5.0 

180304 

JMA-860 

% 

46. /60. 

180304 

JMA-850 

* 

44. /50. 

*  no  illumination  time  given  for  predicting  ERP  level 


The  second  validation  method  involves  comparing  ERP  data  from  ships 
that  make  more  than  one  voyage  through  the  Bay.  Figs.  29  and  30  shows  the 
ERP  variation  for  these  ships.  Approximately  80%  of  the  time,  the  ERP 
calculations  are  within  +/-  1  dB  for  the  repeat  data.  This  is  reasonable,  since 
various  weather  factors  such  as  rain  attenuation,  wave  height,  ducting  and  ship 
wakes  all  influence  the  data  and  were  taken  into  account  only  if  they 
significantly  corrupted  the  ERP  interference  pattern.  There  were  a  few  cases 
where  the  deviations  were  as  high  as  4  dB,  but  these  are  believed  to  be  the 
result  of  situations  with  strong  ducting  conditions  or  anomalous  conditions  such 
as  the  signal  hitting  a  strong  ship  wake  at  the  curve  peak  where  the  ERP  value 
is  taken.  There  is  also  the  possibility  that  the  radar  was  refitted  with  a 
different  transmitter  between  voyages.  The  range  of  ERP  data  was  also 
compared  to  compiled  data  on  ranges  of  ship  radars  and  was  found  to  correlate 
well  (compare  Figures  25  and  31).  The  CBD  data  mainly  includes  merchant 
shipping,  so  does  not  have  as  many  radars  at  lower  ERP  levels  as  the  published 
listing. 
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ERP  Level  (dB) 


Fig.  31  X-band  radars  from  published  listing  (number  of  radars  indicates 

specific  model  types) 


A  more  accurate  prediction  program  could  be  developed  to  analyze  the 
data.  The  current  program  only  uses  the  first  peak  of  the  data  at  the  closest 
distance  to  determine  the  ERP.  This  could  be  expanded  to  perform  a  best  fit 
using  all  the  data.  Also,  weather  conditions  could  be  incorporated  into  the 
model.  Hourly  weather  data  is  available  from  various  stations  close  to  the 
collection  site  (notably  BWI  airport  and  Pax  River  Naval  Air  Station).  This  could 
be  used  to  include  wave  height  data  and  rain  attenuation  into  the  model.  There 
would  be  no  practical  method  of  incorporating  ship  wake  information  or  ducting 
into  the  model  however.  As  seen  from  the  computer  modeling,  though,  ducting 
mainly  influences,  the  data  beyond  10  nm. 


VIII.  SUMMARY 


An  overview  of  sea  surface  multipath  effects  on  accurate  determination  of 
ERP  data  from  ship  radars  was  presented.  It  was  found  that  a  simple  model  of 
multipath  interference  effects  could  be  used  to  quantify  ship  radar  ERPs  over 
short  (5-20  nm)  ranges  to  within  +/-  1  dB  accuracy.  Incorporation  of  weather 
effects  such  as  rain  attenuation  and  wave  height  are  expected  to  provide  some 
improvement  to  the  accuracy,  however,  other  effects  such  as  ship  wakes  and 
ducting  contribute  errors  and  would  be  very  difficult  to  quantify  in  an 
automated  system. 

Several  computer  programs  were  used  to  analyze  the  radar  amplitude 
data.  It  was  found  that  the  EREPS  code  was  only  good  for  coarse  predictions 
because  the  refractivity  profile  is  assumed  to  be  linear.  The  RPO  code  was 
found  to  be  very  useful  for  ground  links.  Its  predictions  for  these  links  are 
thought  to  be  limited  by  the  amount  of  refractivity  data  available.  In  general, 
refractivity  data  with  the  kind  of  resolution  needed  to  accurately  characterize  rf 
propagation  is  difficult  to  obtain  at  even  at  one  location,  much  less  over  the 
entire  length  of  a  link.  In  our  case,  using  one  refractivity  profile  for  the  entire  9 
nm  over-water  link  did  not  provide  exact  predictions  for  link  performance, 
although  relative  performance  estimates  could  be  obtained  and  are  valuable  for 
analyzing  the  link  data. 
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Appendix  A. 
Overview  of  ducting 


A  duct  is  a  layer  of  air  with  different  properties  than  its  surroundings 
which  acts  as  a  waveguide  to  trap  electromagnetic  energy.  Ducting  can  occur 
with  various  meteorological  conditions. 

The  base  of  a  duct  can  be  at  the  sea  surface  or  elevated  above  it.  Ducting 
effects  are  more  prominent  at  longer  ranges,  and  allow  reception  of  RF  signals  at 
much  longer  ranges  than  would  normally  occur.  Fig.  A-l  shows  a  schematic  of 
ducting. 


Ducting  can  cause  deep  fades  in  links.  It  is  highly  variable  and  can  be 
time-of-day  and/or  seasonally  dependent.  It  is  usually  worse  in  the  summer 
months.  Ducting  can  also  cause  location  errors.  Figure  A-2  shows  the  effect  of 
an  evaporation  duct  on  an  over-water  link. 


3  6 
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Fig.  A-2  Propagation  loss  versus  range  for  over-water  link  at  9400  MHz  with 
and  without  evaporation  duct  (EREPS  code,  Wind  speed  =  10  knots,  Transmitter 
height  =  60  ft,  Receiver  height  =180  ft,  no  ducting,  omni  transmitter) 

Low  altitude  propagation  over  water  can  be  influenced  by  surface-based 
ducts  and  evaporation  ducts  [A-l,A-2].  Surface-based  ducts  generally  occur  at 
higher  heights  than  evaporation  ducts.  They  are  caused  by  temperature  and 
humidity  inversions  when  a  warm  dry  air  mass  moves  over  a  cool,  moist  air 
mass.  Surface-based  ducts  are  much  less  common  than  evaporation  ducts, 
though  they  generally  have  larger  thickness. 

Evaporation  ducts  are  a  very  common  occurrence  over  ocean  surfaces. 
They  are  primarily  caused  by  the  rapid  decrease  in  humidity  with  increasing 
altitude  just  above  the  ocean's  surface.  This  humidity  gradient  is  always 
present  over  the  ocean.  The  altitude,  though,  where  the  humidity  achieves  the 
nominal  ambient  value  can  vary  greatly.  This  altitude  is  called  the  duct  height. 
Table  A-I  lists  the  average  duct  heights  for  various  global  regions  [A-l]. 
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Table  A-I  Average  duct  heights  for  various  global  regions  [A-l] 


Area  descriptor 

Evap. 

duct 

heisht.  m 

SFC-based 

duct 

height,  m 

Occurrence 
SFC  duct 

% 

Northern  Atlantic 

5.3 

42 

1.3 

Eastern  Atlantic 

7.4 

64 

2.8 

Canadian  Atlantic 

5.8 

86 

4.1 

Western  Atlantic 

14.1 

118 

9.8 

Mediterranean 

11.8 

125 

13.4 

Persian  Gulf 

14.7 

202 

45.5 

Indian  Ocean 

15.9 

110 

13.4 

Tropics 

15.9 

99 

13.6 

Northern  Pacific 

7.8 

74 

6.2 

Worldwide  average 

13.1 

oo 

U3 

8.0 

For  radar  systems,  ducting  can  cause  a  large  increase  in  the  received 
signal  for  targets  within  the  duct  and  extend  the  radar  detection  range  [A-3, 
A-4].  Radar  holes  can  also  occur  where  returns  would  be  expected  in  non¬ 
ducting  conditions.  There  will  also  be  a  strong  increase  in  sea  clutter. 

In  general,  ducting  effects  occur  at  ranges  longer  than  those  necessary  for 
the  ERP  calculations  in  this  report.  Occasionally,  ducting  will  cause  amplitude 
variations  from  the  free  space  values  at  the  longer  ranges  (10-20  nm),  but  these 
can  be  easily  identified  because  of  the  amplitude  distortion  of  the  multipath 
signal  when  viewed  over  a  range  of  distances. 
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Appendix  B. 
Navigation  Radars 


The  ideal  beam  pattern  for  navigation  radars  is  a  fan  beam,  narrow  in  the 
horizontal  plane  and  wide  in  the  vertical  plane.  A  narrow  beam  is  necessary  to 
produce  high  gain  and  sharp  resolution  as  the  beam  scans  the  horizon  [B-l]. 

The  vertical  pattern  is  broader  so  movements  of  the  ship  do  not  cause  the  beam 
to  miss  the  region  of  the  surface  out  to  the  horizon.  Early  navigation  radars 
used  parabolic  reflector  antennas  to  produce  this  pattern.  For  the  last  20-30 
years,  however,  center  or  end-fed  slotted  waveguide  antennas  have  been 
preferred.  They  have  a  simple  design,  are  reliable,  and  are  easy  to  maintain. 
They  usually  have  a  radome  covered  end  or  are  completely  housed  in  a  radome. 
Longer  length  antennas  have  higher  gain.  The  azimuthal  (horizontal)  3  dB 
beamwidth  of  a  typical  antenna  is  2°.  The  antennas  are  almost  all  horizontally 
polarized.  Horizontal  polarization  does  not  interact  with  the  water  surface  as 
severely  as  vertical  polarization,  so  the  beam  can  scan  the  horizon  more 
effectively. 

Navigation  radar  frequencies  fall  in  the  3  cm  or  10  cm  bands  (X  or 
S  band).  Typical  frequencies  are  9350-9425  MHz  and  2950-3100  MHz. 

The  radar  heights  for  the  ship  transmitters  is  approximately  8-10  m 
(26-33  ft.)  above  sea  level  for  tugs  and  small  boats,  15-20  m  (49-66  ft.)  for 
large  boats  1000-5000  tons,  and  25-30  m  (82-98  ft.)  for  large  boats  5000  tons 
and  above. 
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Appendix  C. 

Calculation  of  Effective  Earth  Radius,  Refractivity,  and 
Altitude  from  Barometric  Readings 


The  effective  earth  radius  can  be  determined  from  [C-l] 


i  = _ _ _ 

,  a'  dn 
n  dh 


where 

a’=  local  radius  of  the  earth  =  sea  level  radius  +  height  above  sea  level 
a  =  sea  level  radius  of  the  earth  =  3963  miles  =  6380  km.  (1.853  km/nm) 
n  =  refractive  index 

dn 

—  =  vertical  gradient  of  the  refractive  index 
dh 

dn!  dh  is  normally  negative  and  on  the  order  of  -40  x  10'6/km. 


N  is  the  refractivity 
N  =  (n-\)\06 

N  =  77.6-  +  3.73  *105-^ 

T  T2 

where  P  is  the  total  atmospheric  pressure  in  millibars, 
e  is  the  water-vapor  pressure  in  millibars 
T  is  the  absolute  temperature  in  degrees  Kelvin. 

The  water-vapor  pressure  is  calculated  as  the  saturation-vapor  pressure 
at  temperature  T,  multiplied  by  the  relative  humidity  (expressed  as  a  decimal). 
[Table  C-l] 

Modified  Refractivity,  M  =  N  +  0.157  h  for  altitude  h  in  meters 
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Table  C-I  Saturation  Vapor  Pressure  over  Water  and  Ice,  in 
millibars  [C-2] 


Temperature,  deg.  C 

Units 

TENS 

0 

1 

2 

3 

4 

5 

6 

7 

8 

9 

40 

73.777 

77.802 

82.015 

86.423 

91.034 

95.855 

100.89 

106.16 

111.66 

117.40 

30 

42.43 

44.927 

47.551 

50.307 

53.2 

56.236 

59.422 

62.762 

66.264 

69.934 

20 

23.373 

24.861 

26.43 

28.086 

29.831 

31.671 

33.608 

35.649 

37.796 

40.055 

10 

12.272 

13.119 

14.017 

14.969 

15.977 

17.044 

18.173 

19.367 

20.63 

21.964 

+0 

6.1078 

6.5662 

7.0547 

7.5753 

8.1294 

8.7192 

9.3465 

10.013 

10.722 

11.474 

-0 

6.1078 

5.623 

5.173 

4.757 

4.372 

4.015 

3.685 

3.379 

3.097 

2.837 

Altitude  can  be  determined  by  knowing  the  difference  in  pressure 
between  two  heights  as  well  as  the  temperature  of  the  air  between  them  [C-2]. 


h2-hi 


RT*  In  p1  —  In  p2 

m  9.8xl04 


where 

=  height  at  position  2,  meters 
hx  =  height  at  position  1,  meters 
p2  —  pressure  at  position  2,  in.  Hg 
px  =  pressure  at  position  1,  in.  Hg 

R=  universal  gas  constant  =  8.3143  x  10^  ergs  K"1  mol'l 

T*  =  temperature  between  two  heights,  deg.  K 
m  =  gram-molecular  weight  of  air  ~  28.9 

Information  on  atmospheric  refractivity  and  archived  refractivity  data  is 
available  from  several  sources  [C-3  -  C-7]. 
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Appendix  D. 

Wind  Speed  versus  Wave  Height: 

The  EREPS  program  uses  wind  speed  to  approximate  wave  height.  The 
wave  height,  Havg,  is  a  function  of  wind  speed,  Ws,  in  meters/sec  according  to 
the  following  equation: 


H 


avg  - 


w< 


8.67 


2.5 


Note:  1  knot  is  1  nautical  mile  per  hour  (1.15  mph).  One  nautical  mile  is  1852  m 
or  6076  feet.  One  knot  is  equal  to  0.514  meters/sec. 


Wind 

Wind 

Wave 

Wave 

speed 

(m/s) 

height 

(m) 

height 

(ft) 

5  knots 

2.57 

0.048 

0.16 

10  knots 

5.14 

0.27 

0.86 

15  knots 

7.71 

0.75 

2.46 

20  knots 

10.28 

1.53 

5.02 
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APPENDIX  E. 
EREPS  Modeling  Results 


Fig.  E-l  Propagation  Loss  versus  range  at  9400  MHz  (EREPS  code,  Wind  speed  10 
knots,  Transmitter  height  =  60  ft,  Receiver  height  =  180  ft,  no  ducting,  omni 

transmitter) 
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Range,  nautical  miles 

Fig.  E-2  Propagation  Loss  versus  range  at  3050  MHz  (EREPS  code,  Wind  speed  10 
knots.  Transmitter  height  =  60  ft,  Receiver  height  =180  ft,  no  ducting,  omni 

transmitter) 


Range,  nautical  miles 

Fig.  E-3  Propagation  Loss  versus  range  at  9400  MHz  with  varying  receiver  heights 
(EREPS  code,  Wind  speed  =10  knots,  Transmitter  height  =  60  ft,  no  ducting,  omni 

transmitter) 


Fig.  E-4  Propagation  Loss  versus  range  at  9400  MHz  with  varying  receiver  heights 
(EREPS  code,  Wind  speed  =  10  knots,  Transmitter  height  =  60  ft,  no  ducting,  omni 

transmitter) 


Range,  nautical  miles 

Fig.  E-5  Propagation  Loss  versus  range  at  9400  MHz  with  varying  receiver  heights 
(EREPS  code,  Wind  speed  =  10  knots,  Transmitter  height  =  60  ft,  no  ducting,  omni 

transmitter) 
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5  20 

Range,  nautical  miles 

Fig.  E-6  EREPs  Propagation  Loss,  9400  MHz  with  double  receiver  (150  and  180  ft) 
(Wind  speed  =10  knots,  Transmitter  ht.  =  60  ft,  no  ducting,  omni  transmitter) 
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Range,  nautical  miles 

Fig.  E-7  Propagation  Loss  versus  range  at  9400  MHz  with  varying  receiver  heights 
(150  and  180  ft)  (EREPS  code,  Wind  speed  =10  knots.  Transmitter  height  =  80  ft,  no 

ducting,  omni  transmitter) 


Range,  nautical  miles 

Fig.  E-8  Propagation  Loss  versus  range  at  9400  MHz  with  combined  receiver  (150 
and  180  ft)  (EREPS  code,  Wind  speed  =  10  knots,  Transmitter  height  =  80  ft,  no  ducting, 

omni  transmitter) 


Range,  nautical  miles 

Fig.  E-9  Propagation  Loss  versus  range  at  3050  MHz  with  varying  receiver  heights 
(150  and  180  ft)  (EREPS  code.  Wind  speed  =  10  knots,  Transmitter  height  =  60  ft,  no 

ducting,  omni  transmitter) 
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APPENDIX  F 

Codes  for  Multipath  Modeling* 

1.  Accelerated-Intermediate-Region  Propagation  Model 
Dave  Kenney  (JSC),  (410)  573-7371 

2.  CORAFIN  (France) 

Eric  Mandine  (Minister  de  la  Defense  University  de  Toulon) 

Phone:  94  11  40  47  (poste  29  047) 

3.  DIF-CERT  (France) 

Douchin  (?) 

4.  EEMS  1.0  (Hybrid  model  that  supersedes  PCPEM,  Ver.  1.0)  -  UK 

Mireille  Levy  (Radio  Comms.  Research  Unit)  +44  (0)  1235  446522 
Kenneth  Craig  ("  ’’  "  ")  +44  (0)  1235  445134 

5.  EREPS/IREPS  (available  over  the  internet) 

Clause  Hatten  (NRaD),  (619)  553-1427 

6.  GELTI 

Ray  Luebbers 

7.  HAP  (Israel):  Sherman  Marcus  (RAFAEL) 

8.  IDFG:  Sherman  Marcus  (RAFAEL) 

9.  MLAYER  (Multi-Layer  Waveguide) 

Ken  D.  Anderson  (NRad),  (619)  553-1420 

10.  PCPEM  (UK) 

1 1 .  RPO  (Radio  Physical  Optics)  (available  over  the  internet) 

Herb  V.  Hitney  (NRad)  (619)  553-1428 

12.  TEMPER  (Tropospheric  Electromagnetic  Parabolic  Equation  Routine) 

Dan  Dockery  (JHU/APL)  (301)  953-5461 

13.  Tirem  (Terrain  Integrated  Rough  Earth  Model) 

David  W.  Eppink  (IITRI/JSC)  (410)  573-7599 

14.  TPEM  (Terrain  Parabolic  Equation  Model)  (available  over  the  internet) 

Amalia  Barrios  (NRad)  (619)  553-1429 

15.  VTRPE  (Variable  Terrain  Parabolic  Equation  Model) 

Frank  J.  Ryan  (NRad)  (619)  553-3099 

*A  majority  of  the  codes  on  this  list  were  obtained  from  the  Electromagnetic  Workshop 
at  the  Applied  Physics  Lab  in  July  1995 
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Appendix  G. 

ERP  Determination  from  Ship  Data: 
Sample  Plots 


Fig.  G-l  S-band,  3049  MHz,  Radid  342157,  out-bound,  12/8/97,  1700,  Transmitter 
height  121  feet,  receiver  height,  116  ft.  Free  space  loss  is  approximately  5  dB  below 
the  peak  signal  (varies  slightly  with  sea  state).  The  offset  for  the  free  space  loss 
corresponds  to  the  ERP  of  the  generating  signal,  in  this  case  105.5  dBm. 
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Distance,  nm 

Fig.  G-2  S-band,  3051  MHz,  Radid  309325,  in-bound,  11/05/97,  2200,  Transmitter 
height  139  feet,  receiver  height,  116  ft.  Free  space  loss  is  approximately  5  dB  below 
the  peak  signal  (varies  slightly  with  sea  state).  ERP  =  103.  dBm. 


Distance,  nm 

Fig.  G-3  S-band,  3053  MHz,  Radid  311496,  out-bound,  11/07/97,  1100,  Transmitter 
height  111  feet,  receiver  height,  116  ft.  Free  space  loss  is  approximately  5  dB  below 
the  peak  signal  (varies  slightly  with  sea  state).  ERP  =  103.  dBm. 


Distance,  nm 

Fig.  G-4  S-band,  3058  MHz,  Radid  154860,  in-bound,  6/3/97,  1200,  Transmitter 
height  115  feet,  receiver  height,  116  ft.  Free  space  loss  is  approximately  5  dB  below 
the  peak  signal  (varies  slightly  with  sea  state).  ERP  =  103.  dBm. 


Distance,  nm 

Fig.  G-5  X-band,  9363  MHz,  Radid  342157,  out-bound,  12/8/97,  1700,  Transmitter 
height  82  feet,  receiver  height,  116  ft.  Free  space  loss  is  approximately  5  dB  below  the 
peak  signal  (varies  slightly  with  sea  state).  ERP  =  103.1  dBm. 
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Signal  amplitude,  dB 


APPENDIX  H. 

Upper  Air  Soundings 

Figures  H-l  and  H-2  show  sounding  data  from  6/17/96.  This  data  was  used  in 
the  RPO  code  to  produce  the  propagation  loss  diagrams  shown  in  Figures  H-3  and  H-4. 
Figures  H-5  and  H-6  show  propagation  loss  diagrams  for  two  RPO  runs  with  standard 
linear  atmospheric  refractivity  profiles. 

Figures  H-7  through  H-13  show  sounding  data  from  7/16/96.  This  data  was  used 
in  the  RPO  code  to  produce  the  propagation  loss  diagrams  shown  in  Figures  H-l 4 
through  H-17.  Figure  H-l 8  shows  a  propagation  loss  diagram  for  one  RPO  run  with  a 
standard  linear  atmospheric  refractivity  profiles. 

Figures  H-19  and  H-20  show  the  Tilghman  Island  Radar  Signal  Level  versus  Time 
for  the  period  during  7/16/96  when  the  soundings  data  was  taken.  Tide  Height  is  also 
shown. 


Height,  m 

Fig.  H-l  Upper  air  sounding,  6/17/96,  Set  #1,  10:30  am,  Weather  Conditions:  sunny, 
hot,  winds  6  knots  at  10:30  am,  10-13  knots  at  1:30  pm 
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■  Temperature,  deg.  F 

■  Rel.  Humidity  % 


Fig.  H-2  Upper  air  sounding,  6/17/96,  Set  #2,  1:30  pm,  Weather  Conditions:  sunny 
hot,  winds  6  knots  at  10:30  am,  10-13  knots  at  1:30  pm 


Distance,  km 


110  115  120  125  130  135  140  145  150 
Propagation  Loss,  dB 

Fig.  H-3  RPO  run  using  sounding  data  from  6/17/96,  Set  #1,  10:30  am,  tide  ht.=l  ft. 


0  5  10  15  20 

Distance,  km 


110  115  120  125  130  135  140  145  150 
Propagation  Loss,  dB 

Fig.  H-4  RPO  run  using  sounding  data  from  6/17/96,  Set  #2,  1:30  pm,  tide  ht.=  l  ft. 


0  5  10  15  20 

Distance,  km 


110  115  120  125  130  135  140  145  150 
Propagation  Loss,  dB 

Fig.  H-5  RPO  run,  stdatm,  tide  height  =  1  ft.,  0  350,  1000.  468 
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Height, 


0  5  10  15  20 

Distance,  km 


110  115  120  125  130  135  140  145  150 
Propagation  Loss 


Fig.  H-6  RPO  run,  stdatm,  tide  height  =  1  ft.,  0  350,  1000.  420  note:  lines  get  closer 

together 


r* 
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Relative  Humidity  %  or  Temperature,  deg.  F 


110 


Fig.  H-7  Upper  Air  Sounding,  7/16/96  Set  #1,  10:20  am 


360  380  400  420  440 

Modified  Refractivity,  M-units 

Fig.  H-8  Upper  Air  Sounding,  7/16/96  Set  #1,  10:20  am 
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Relative  Humidity  %  and  Temp.  deg.  F 


Temperature,  deg.  F,  and  Relative  Humidity  % 
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Fig.  H-13  Upper  Air  Sounding,  7/16/96,  Comparing  Modified  Refractivity  Profile  of 
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110  115  120  125  130  135  140  145  150 
Propagation  Loss,  dB 

Fig.  H-14  RPO  run  using  7/16/96  Sounding  Data,  Set  #1,  10:20  am,  rpo  code,  tide 

height  0  ft. 
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Distance,  km 
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Fig.  H-15  RPO  run  using  7/16/96  Sounding  Data,  Set  #2,  12:20  pm,  rpo  code  run,  tide 
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Fig.  H-16 


^  I  * 

110  115  120  125  130  135  140  145  150 
Propagation  Loss,  dB 


RPO  run  using  7/16/96  Sounding  Data,  Set  #3,  2:00  pm,  rpo 

height  0 


code  run,  tide 
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0  2  4  6  8  10  12  14  16  18 

Distance,  km 


Fig.  H-17 


& 


•P 

b> 


•H 

CD 


lljpi 

110  115  120  125  130  135  140  145  150 
Propagation  Loss,  dB 


RPO  run  using  7/16/96  Sounding  Data,  Set  #3,  2:00  pm,  rpo  code  run,  tide 
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Propagation  Loss,  dB 

Fig.  H-18  RPO  code  run  with  standard  atmospheric  case,  (2  data  pts.)  tide  height=0 
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Fig.  H-19  Tilghman  Island  Radar  Signal  Level  and  Tide  Height,  7/16 


16  16.2  16.4  16.6  16.6  17 

Time,  days 

Fig.  H-20  Tilghman  Island  Radar  Signal  Level  and  Tide  Height,  7/16/96 
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APPENDIX  I. 
Weather  Data  Sources 


Various  sources  are  available  over  the  internet  for  obtaining  weather  data.  This 
data  is  usually  presented  on  a  real  time  basis,  so  archiving  is  necessary  to  use  old  data. 
Tide  data,  radiosonde  data,  weather  station  data,  satellite  visible  and  infrared  (IR) 
pictures  (GOES  satellite),  and  wave  height  and  water  vapor  data  (TOPEX  satellite)  are 
available. 

Temperature,  wind  speed,  wind  direction,  humidity,  water  temperature,  dew 
point,  and  barometric  pressure  are  weather  data  taken  daily  at  CBD.  Weather 
information  about  the  Bay  can  also  be  accessed  over  the  telephone  from  the  weather 
bureau.  By  calling  703-260-0505,  a  forecast  is  given  for  the  Chesapeake  Bay  area  that 
includes  temperature,  wind  speed  and  direction,  and  wave  height. 

Moon  phase  information  is  located  in  the  local  paper  (Washington  Times  or 
Washington  Post).  They  also  give  high  tide  times. 

Tide  data  can  be  accessed  by  telneting  to  WLNET2.NOS.NOAA.GOV.  The  word 
GUEST  is  used  as  a  username.  Archived  six-minute  water  level  data  can  then  be 
accessed  from  various  stations.  Local  stations  are  at  Annapolis  and  Solomon's  Island. 
Tide  data  for  CBD  can  then  be  approximated  by  subtracting  one  hour  from  the 
Annapolis  data.  (Tide  map  shows  1  hour  difference). 

The  Northeast  Regional  Climate  Center  maintains  archived  weather  data  for 
numerous  sites  including  the  Baltimore- Washington  International  airport  (BWI)  and 
the  Patuxant  Naval  Air  Station.  Archived  data  is  obtained  with  a  small  fee.  Daily  data 
is  also  available  and  can  be  archived  independently.  The  internet  address  is 
http://met-www.cit.cornell.edu/nrcc_home.html. 

Topex  satellite  data  including  wave  height  and  water  vapor  content  can  be 
accessed  with  the  internet  address: 

http://podaac-www.jpl.nasa.gov/topex/ 


Michigan  State  University  has  current  weather  maps  and  movies  available  for 
downloading.  They  make  composite  weather  maps  from  GOES  satellite  data.  Updates 
are  every  3  and  6  hours.  Of  interest  to  this  project  is  the  Worldwide  Montage  and 
Worldwide  IR  Composite  pictures.  The  internet  address  is: 
http://rs560.cl.msu.edu/weather/ 

Another  good  site  for  weather  information  is  the  WeatherNet  WeatherSites  page 
on  the  internet.  It  gives  a  list  of  lots  of  weather  sites. 

http://cirrus.sprl.umich.edu/wxnet/servers.html 

Radiosonde  data  is  presented  at  various  sites.  Use  search  engine  under 
'radiosonde'. 
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Weather  archiving  has  been  ongoing  since  January  1996.  This  data  can  be  used  in 
the  analysis  of  open-ocean  tracks.  The  archiving  of  data  off  the  internet  was 
automated.  It  includes  whole  world  satellite  visible  and  IR  imagery  as  well  as  Doppler 
radar  information  where  available.  Hourly  weather  station  data  is  also  archived  for 
various  sites  in  the  DELMARVA  area.  Figure  1-1  shows  GOES  8  data.  Figure  1-2  shows 
a  composite  of  satellite  and  Doppler  radar  images. 
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Moderate  storm: 
overcast  and  moderate 
rainfall 


Large  storm:  at  center  of 
storm,  heavy  rainfall  and 
downpours  likely 

Fig.  1-1  GOES-8  Image 
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